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41 Abstract 

42 

43 To identify host factors relevant for SARS-coronavirus (SARS-CoV) replication, we performed an siRNA 

44 library screen targeting the human kinome. Protein kinases are key regulators of many cellular functions 

45 and the systematic knockdown of their expression should provide a broad perspective on factors and 

46 pathways promoting or antagonizing coronavirus replication. In addition to 40 proteins that promote 

47 SARS-CoV replication, our study identified 90 factors exhibiting an antiviral effect. Pathway analysis 

48 grouped subsets of these factors in specific cellular processes, including the innate immune response and 

49 the metabolism of complex lipids, which thus appear to play a role in SARS-CoV infection. Several 

50 factors were selected for in-depth validation in follow-up experiments. In cells depleted for the |32 subunit 

51 of the coatomer protein complex (COPB2), the strongest proviral hit, we observed reduced SARS-CoV 

52 protein expression and a >2-log reduction in virus yield. Knockdown of the COPB2-related proteins 

53 COPB1 and Golgi-specific brefeldin A-resistance guanine nucleotide exchange factor 1 (GBF1) also 

54 suggested that COPI-coated vesicles and/or the early secretory pathway are important for SARS-CoV 

55 replication. Depletion of the antiviral double-stranded RNA-activated protein kinase (PKR) enhanced 

56 virus replication in the primary screen, and validation experiments confirmed increased SARS-CoV 

57 protein expression and virus production upon PKR depletion. In addition, cyclin dependent kinase 6 

58 (CDK6) was identified as a novel antiviral host factor in SARS-CoV replication. The inventory of pro- 

59 and antiviral host factors and pathways described here substantiates and expands our understanding of 

60 SARS-CoV replication and may contribute to the identification of novel targets for antiviral therapy. 

61 
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62 

63 Importance 

64 

65 Replication of all viruses including SARS-coronavirus (SARS-CoV) depends on and is influenced by 

66 cellular pathways. Although substantial progress has been made in dissecting the coronavirus replicative 

67 cycle, our understanding of the host factors that stimulate (proviral factors) or restrict (antiviral factors) 

68 infection remains far from complete. To study the role of host proteins in SARS-CoV infection, we set out 

69 to systematically identify kinase-regulated processes that influence virus replication. Protein kinases are 

70 key regulators in signal transduction, control a wide variety of cellular processes, and many of them are 

71 targets of approved drugs and other compounds. Our screen identified a variety of hits and will form the 

72 basis for more detailed follow-up studies that should contribute to a better understanding of SARS-CoV 

73 replication and coronavirus-host interactions in general. The identified factors could be interesting targets 

74 for the development of host-directed antiviral therapy to treat infections with SARS-CoV or other 

75 pathogenic coronaviruses. 

76 
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77 Introduction 

78 

79 Positive-stranded RNA (+RNA) viruses interact with the infected host cell at many levels during their 

80 replicative cycle, and thus far numerous host cell proteins that influence vims infection have been 

81 identified (1, 2). These include, for example, host factors recruited by the vims during the various stages 

82 of its replicative cycle and those involved in the host’s defense against vims infection. Such proteins may 

83 constitute interesting targets for the development of novel antiviral strategies, as dmg resistance is less 

84 likely to develop when cellular rather than viral functions are targeted. Antiviral dmg resistance is a 

85 serious problem, in particular when combating RNA vimses, due to their high mutation rate and potential 

86 for rapid adaptation. 

87 Systems biology approaches have been instmmental in advancing our knowledge of the proteins 

88 and cellular pathways that influence +RNA vims infection. For example, systematic functional genomics 

89 screens using small interfering RNA (siRNA) libraries have identified numerous host proteins with a role 

90 in the replication of important human pathogens like West Nile vims (3), Dengue vims (4, 5), human 

91 immunodeficiency vims 1 (6), hepatitis C vims (7-12), and influenza vims (8, 13, 14). For coronavimses a 

92 number of relevant host proteins have previously been described ((15-17), and reviewed in (2, 18)), but 

93 the use of siRNA screens to systematically identify such factors has not been reported thus far. 

94 Coronavimses, and some other members of the order Nidovirales (19), have the largest RNA 

95 genomes known to date (25-34 kb (20)) and the complexity of their molecular biology clearly 

96 distinguishes them from other +RNA vims groups. Although infection with most established human 

97 coronavimses is associated with relatively mild respiratory symptoms (21, 22), the 2003 outbreak of 

98 severe acute respiratory syndrome (SARS) highlighted the potential of zoonotic coronavimses to cause 

99 lethal disease in humans. The emergence of SARS-coronavims (SARS-CoV), which likely originated 

100 from bats, initiated an outbreak that affected about 8,000 humans, with a mortality rate of approximately 

101 10% (23). Strikingly, a similar outbreak of coronavims-induced severe respiratory disease has been 
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102 developing in a number of Arab countries since April 2012, with ~420 of the >1100 confirmed cases 

103 having a fatal outcome by April 2015 (http://www.who.int/). The causative agent, Middle East Respiratory 

104 Syndrome-coronavirus (MERS-CoV), was identified as a previously unknown member of the 

105 betacoronavirus subgroup 2c (24, 25). These recent developments stress the importance of developing 

106 antiviral approaches to combat coronavirus infections and highlight the relevance of the systematic 

107 dissection of coronavirus-host interactions. 

108 SARS-CoV RNA synthesis, like that of many +RNA viruses (26), takes place at virus-induced 

109 membrane structures (27, 28), which in this case comprise a reticulovesicular network (RVN) of modified 

110 endoplasmic reticulum ((28) and reviewed in (29)). The viral replication and transcription complexes 

111 (RTCs) are associated with this RVN, which is thought to create a suitable microenvironment for RNA 

112 synthesis and possibly also provides protection against cellular antiviral activities. The biogenesis of the 

113 RVN, and the functional details of the RTC, in particular the role of cellular factors and pathways, are far 

114 from understood. 

115 Previous studies addressed coronavirus-induced immune responses, as well as a number of 

116 specific interactions between coronaviruses and the antiviral immune response (reviewed in (2)). Several 

117 immune evasion mechanisms were attributed to protein functions that are either conserved across Co Vs or 

118 specific for certain CoV lineages. Proteins such as non-structural protein 1 (nspl; (30), the nsp3 papain- 

119 like proteinase (31), the nspl6 2’-0-methyltransferase (32), the nucleocapsid (N) protein (33), and the 

120 products of SARS-CoV ORFs 3b and 6 (34-37) have been reported to interfere with interferon (IFN) 

121 induction and/or signalling. In addition, the SARS-CoV E protein has been shown to manipulate the 

122 cellular stress response in cell culture, including the unfolded protein response and apoptosis (38). 

123 To gain more insight into the role of host factors in the SARS-CoV replicative cycle, we set out to 

124 systematically identify kinase-regulated cellular processes that influence virus replication. Protein kinases 

125 are key regulators in signal transduction and control a wide variety of cellular processes. Thus, assessing 

126 their relevance for virus replication can provide a broad perspective on factors and pathways relevant for 
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127 SARS-CoV replication, as illustrated by previous studies identifying cellular kinases as host factors 

128 influencing various stages of the replicative cycle of other +RNA viruses (5, 10, 11, 39, 40). 

129 In this study, we have screened an siRNA library that targets the cellular kinome (779 genes) and 

130 identified 40 proviral and 90 antiviral factors whose depletion significantly reduced or enhanced SARS- 

131 CoV replication, respectively. Pathway analysis grouped several subsets of hits in specific cellular 

132 pathways, suggesting that these play an important role in the SARS-CoV-infected cell. Two strong hits 

133 from the siRNA screen, the proviral (32 subunit of the coatomer complex (COPB2) and the antiviral 

134 double-stranded RNA-activated protein kinase (PKR), were selected for independent validation and 

135 follow-up analysis, which confirmed their importance for SARS-CoV replication. In addition, several 

136 other hits from the primary screen were evaluated, and the relevance of the antiviral factor CDK6 and the 

137 proviral factor PRKCi could be confirmed. Our data offer a glimpse into the complex interplay between 

138 SARS-CoV and its host cell, and provide a basis for in-depth studies that will enhance our understanding 

139 of coronavirus replication and coronavirus-host interactions. 

140 

141 Materials and methods 

142 

143 Cell culture, compound, viruses, and virus titration - 293/ACE2 (41) and Vero E6 cells were cultured 

144 as described previously (42). Although 293/ACE2 cells have been described as a human 293 cell-derived 

145 cell line (41), our recent work established that these cells actually must have originated from a non-human 

146 primate species that is closely related to the rhesus monkeys Macaca mulatto and Papio Anubis (43). Cells 

147 were infected with SARS-CoV strain Frankfurt-1 (44) or GFP-expressing recombinant SARS-CoV 

148 (Urbani strain) (45) as described previously (42). Sodium aurothiomalate (ATM; Sigma cat. nr. 157201) 

149 was dissolved in PBS and stored as 100 mM stock at -20°C.Virus titrations were performed essentially as 

150 described before (46). All work with infectious wild-type (wt) SARS-CoV and SARS-CoV-GFP was 

151 performed inside biosafety cabinets in a biosafety level 3 facility at Leiden University Medical Center. 
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152 

153 siRNA library and transfection reagents - The ON-TARGETplus SMARTpool Protein Kinases siRNA 

154 Library that targets the itlRNAs of 779 genes, comprising the complete human kinome and some 

155 additional targets, was obtained from Dharmacon. Each individual siRNA SMARTpool consisted of four 

156 siRNAs targeting the same gene. A non-targeting (scrambled) siRNA (cat. nr. D-001810-10; Dharmacon) 

157 served as a negative control and a GAPDEI-targeting siRNA (cat. nr. D-001830-10; Dharmacon) was used 

158 to routinely monitor transfection and knockdown efficiency. Stock solutions (2 pM) of siRNA 

159 SMARTpools were prepared by dissolving 0.5 nmol of an siRNA SMARTpool in 250 pi of lx siRNA 

160 buffer (Dharmacon), according to the manufacturer’s instructions. Using a 96-well pipettor (Rainin 

161 Liquidator 96), the contents of the siRNA library master plates was aliquoted into volumes appropriate for 

162 individual screening experiments. The resulting sets of ten deep-well 96-well library plates (Greiner Bio- 

163 One) were stored at -80°C until further use. 

164 

165 siRNA library screening and validation - In each siRNA screen, 293/ACE2 cells in 96-well plates 

166 containing ~10 4 cells per well were transfected with a 100-pl mixture containing 100 nM siRNA, 0.2 pg 

167 DharmaFECTl (Dharmacon), OptiMEM (Invitrogen), and antibiotic-free cell culture medium, 

168 supplemented with 8% fetal calf serum (FCS) and 2.5 mM L-Glutamine, according to Dharmacon's 

169 instructions. Transfection mixes were prepared in the ten deep-well 96-well plates that together contained 

170 the complete library of 779 siRNA SMARTpools (see above). Using the contents of these library plates, 

171 we transfected 293/ACE2 cells in black (3 wells per target) and transparent 96-well plates (3 wells per 

172 target). For a schematic representation of the experimental set-up, see Fig. 2. Transfection of individual 

173 siRNAs (ON-TARGETplus siRNAs; Dharmacon) targeting CDK6 (cat. nr. LU-003240-00), MAP2K1 

174 (cat. nr. LU-003571-00), MAP2K3, (cat. nr. LU-003509-00), PKR (cat. nr. LU-003527-00), or siRNA 

175 SMARTpools targeting COPB1 (cat. nr. L-017940-01) and GBF1 (cat. nr. L-019783-00) was performed 

176 as described previously (42). Twenty-four hours post transfection (p.t.), the medium was replaced, and 
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177 cells were incubated for another 24 h at 37°C. At 48 h p.t., cells were infected with SARS-CoV-GFP at an 

178 MOI of 10, and 24 h later they were fixed with 3% paraformaldehyde (PFA) in PBS. GFP expression was 

179 quantified by measuring fluorescence in a 96-well plate reader (Berthold Mithras LB 940), using 

180 excitation and emission wavelengths of 485 and 535 nm, respectively. The fluorescence in wells 

181 containing mock-infected cells was used to correct for background signal. 

182 

183 GAPDH and cell viability assays - At 48 h p.t., GAPDFI enzyme activity in lysates of siRNA-transfected 

184 cells was measured using the KDalert™ GAPDFI Assay Kit (Ambion) according to the manufacturer's 

185 instructions. Possible cytotoxic effects of siRNA transfection were analyzed (in triplicate) at 48 h p.t., 

186 using the CellTiter 96® AQ ueous Non-Radioactive Cell Proliferation Assay (Promega). After 90 min, the 

187 reaction was terminated by the addition of 25 pi of 10% SDS and absorbance at 490 nm (A 490 ) was 

188 measured using a 96-well plate reader (Berthold). 

189 

190 Data analysis - Raw data from GFP fluorescence and cell viability measurements were analyzed per 

191 individual screen with the Bioconductor/R package CellHTS2 (47) with minor modifications (see Results 

192 section for details). Average GFP expression (n=3) and cell viability were calculated and normalized to 

193 the signals of scrambled siRNA-transfected (control) cells. A two-sided one-sample Student’s t test was 

194 used on the logi-transformed normalized values to determine the significance (p < 0.05) of the changes in 

195 GFP expression caused by siRNA transfection. The siRNA transfection was considered non-cytotoxic 

196 when the normalized cell viability assay readings (A 490 ) were above 0.85 (p < 0.05). Significance was 

197 determined using a one-sided one-sample Student’s t test on the log 2 -transformed normalized values using 

198 p < 0.85 as the null hypothesis. 

199 

200 Gene silencing using lentivirus-expressed shRNAs - Vectors for expression of short hairpin RNAs 

201 (shRNAs) targeting human COPB2 (cat. nr. TRCN-065114; accession nr. NM 004766) or expression of a 
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202 non-targeting (scrambled) control shRNA (cat. nr. SHC-002) were picked from the MISSION TRC-1 

203 library of shRNA-expressing lentiviruses (Sigma) and lentivirus stocks were prepared according to the 

204 manufacturer’s instructions. Lentivirus particle titers were determined using a p24 ELISA (Zeptometrix) 

205 according to the manufacturer's instructions. Wells (4 cm 2 ) containing 8 x 10 4 293/ACE2 cells were 

206 transduced with shRNA-expressing lentiviruses at an MOI of 3 in culture medium containing 8 pg/ml 

207 polybrene, and after 24 h fresh medium was given. At 72 h p.t., cells were infected with wt SARS-CoV or 

208 SARS-CoV-GFP (MOI 0.01), and depletion of COPB2 was validated by Western blotting. 

209 

210 Protein analysis and antibodies - Total cell lysates were prepared in 4x Laemmli sample buffer (100 

211 mM Tris-HCl, pH 6.8, 40% glycerol, 8% sodium dodecyl sulfate (SDS), 40 mM DTT, 0.04 mg/ml 

212 bromophenol blue), after which samples were heated at 95°C for 15 min. Following SDS-PAGE, proteins 

213 were transferred to Hybond-LFP membranes (GE Healthcare) by semi-dry blotting, and membranes were 

214 blocked with 1% casein in PBS containing 0.1% Tween-20 (PBST). The following antisera against 

215 cellular proteins were used: rabbit anti-PKR (cat. nr. 610764; BD Biosciences), goat anti-COPB2 (sc- 

216 13332; Santa-Cruz), rabbit anti-CDK6 (sc-177; Santa Cruz), rabbit anti-MAP2Kl (710446; Life 

217 Technlogies), rabbit anti-MAP2K3 (sc-961; Santa Cruz), and mouse monoclonal antibodies against (3- 

218 actin (A5316; Sigma) and the transferrin receptor (TfR; cat. nr. 13-6890; Invitrogen). Rabbit antisera 

219 against SARS-CoV nsp8 and N protein (28, 48) were used to analyze viral protein expression. After 

220 overnight incubation with the primary antibody, membranes were probed with biotinylated secondary 

221 antibodies (rabbit anti-goat, swine anti-rabbit, or goat anti-mouse) for 1 h at RT, after which a tertiary 

222 mouse anti-biotin-Cy3 antibody was used to visualize protein bands using a Typhoon 9410 scanner (GE 

223 Healthcare). 

224 

225 Canonical pathway analysis - The Ingenuity Pathway Analysis (IPA™) package was used to assign hits 

226 to canonical cellular pathways. The significance of the association between the dataset and the respective 
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227 pathways was determined in two ways: (i) the number of molecules from the dataset that mapped to a 

228 specific pathway divided by the total number of molecules in that canonical pathway (the higher the 

229 percentage of hits identified in a specific pathway, the higher the likelihood it plays a role in the viral 

230 replicative cycle); (ii) Fisher’s exact test was used to determine the probability that the association 

231 between the genes in the dataset and the canonical pathway is explained by chance alone. 

232 
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233 Results 

234 

235 Development of an siRNA screening protocol for host factors involved in nidovirus replication. A 

236 robust protocol was developed to assess the effect of systematic knockdown of individual host kinases on 

237 the replicative cycle of SARS-CoV (this study) and the distantly related arterivirus equine arteritis virus 

238 (EAV; K. F. Wannee, A. H. de Wilde et al., manuscript in preparation), which was applied to the 

239 screening of a commercial human kinome-directed siRNA library (779 targets). We performed our siRNA 

240 screens in 293/ACE2 cells (41), which express the SARS-CoV receptor angiotensin-converting enzyme 2 

241 and, in contrast to other cell lines tested, were found to be permissive to a combination of siRNA 

242 transfection and infection with either SARS-CoV or EAV. This property facilitated direct comparative 

243 studies between these two distantly related nidoviruses. Unfortunately, after completion of the siRNA 

244 screens, it was discovered that these cells are not of human origin, but have most likely originated from an 

245 Old World monkey closely related to Papio Amtbis and Macaca mulatto (43). Nevertheless, because the 

246 sequence identity between the human genome and that of several Old World monkeys is 94% (49) and 

247 because pools of four siRNAs were used for each target, we believe that the consequences of the 

248 misidentification of this cell line are limited, although the chance of false-negative hits may have been 

249 somewhat increased (49). Infection of 293/ACE2 cells with SARS-CoV-GFP at an MOI of 10 yielded a 

250 robust and readily detectable GFP signal at 24 h p.i. (Fig. 1 A). The GFP signal was stronger at 28 and 30 h 

251 p.i., indicating that it had not yet reached a plateau at 24 h p.i. (Fig. 1A). We therefore chose the latter 

252 time point to fix cells and measure GFP fluorescence, as it should also allow the identification of antiviral 

253 factors whose knockdown would increase reporter gene expression. The 293/ACE2 cells could be 

254 efficiently transfected with siRNAs, as illustrated by a consistent ~75% reduction of GAPDH activity at 

255 48 h p.t. using an siRNA SMARTpool targeting the GAPDH mRNA (Fig. IB; white bars). No change in 

256 cell viability was detected by 48 h p.t. following transfection with either a scrambled siRNA or the 

257 GAPDH-specific siRNA (Fig. IB; grey bars). When these cells were subsequently infected with SARS- 
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258 CoV-GFP (MOI 10), no significant differences in GFP expression were observed at 24 h p.i. compared to 

259 control cells that had not been transfected with siRNAs. This demonstrated that the siRNA transfection 

260 procedure per se did not adversely affect SARS-CoV-GFP replication (Fig. IB; black bars). 

261 

262 Kinome-wide siRNA screens for host factors involved in SARS-CoV replication. A human kinome- 

263 directed siRNA screen was performed to identify host cell kinases that affect SARS-CoV-GFP replication, 

264 according to the experimental set-up outlined in Fig. 2. For each independent siRNA screening 

265 experiment, we used a set of ten 96-well library plates, each containing approximately 80 specific siRNA 

266 SMARTpools and several controls. Transfection mixes (final concentration of 100 nM siRNA) were 

267 prepared in these library plates and their contents was used to transfect - per library plate - 293/ACE2 cells 

268 in three black and three transparent 96-well plates. Forty-eight hours after siRNA transfection, the black 

269 plates were infected with the SARS-CoV-GFP reporter virus (MOI 10), and at 24 h p.i. GFP expression 

270 was measured by fluorometry. At the moment of infection, the transparent plates were used to monitor 

271 (potential) cytotoxic effects of siRNA transfection using a colorimetric cell viability assay. The complete 

272 siRNA screen, i.e. the viability controls (in triplicate for each siRNA SMARTpool) and the quantitation of 

273 SARS-CoV-driven GFP expression (in triplicate), was repeated in three independent experiments. The 

274 data, obtained from a 96-well plate reader, was processed with the Bioconductor/R package CellHTS2 as 

275 described (47). Experimental controls were assigned, and the NPI method (normalized percent of 

276 inhibition) was used to normalize GFP fluorescence values to those of scrambled siRNA-transfected cells, 

277 and to correct for plate-to-plate variation. Subsequently, the GFP data were transformed to a multiplicative 

278 scale (the value obtained using scrambled siRNA-transfected cells was set to 1). Next, the results for each 

279 replicate library screen were summarized and used for further data analysis, including the assignment of 

280 GenelDs to each well. Finally, the data of the three independent library screens were combined and 

281 summarized. 
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282 Host cell kinases were considered to have a proviral effect when their siRNA-mediated 

283 knockdown reduced the GFP signal (negative score values) and kinases were considered antiviral when 

284 the GFP signal increased upon their knockdown (positive score values). Graphical representations of the 

285 hit distribution per plate were visually inspected in order to minimize the chance of false positive or false 

286 negative hits due to major (technical) artifacts (data not shown). 

287 Using scrambled siRNA-transfected control cells as a reference, the knockdown of most cellular 

288 kinases was found to be non-cytotoxic within the time frame of this experiment (Fig. 3A and Dataset SI). 

289 Transfection of siRNAs was considered to be cytotoxic when the viability of cells transfected with a 

290 target-specific siRNA pool was <85% of the viability of control cells transfected with scrambled siRNAs 

291 (Fig. 3A). Using this criterion, 222 out of 779 (28.5%) transfections with the specific siRNA pools 

292 appeared to be toxic to the cells. A minor fraction (50 targets; 6.4%) appeared to be highly detrimental 

293 (normalized viability value below 75%). To prevent false-positive proviral hits due to a general negative 

294 effect on cell viability or cell division, we excluded all targets whose knockdown was associated with 

295 viability measurements below 85%. Such data filtering was not applied for antiviral hits (i.e. hits whose 

296 knockdown enhanced GFP expression) since siRNA-induced cytotoxicity is expected to inhibit virus 

297 replication and should therefore not give rise to false-positive antiviral hits. 

298 

299 Proviral and antiviral proteins and pathways in SARS-CoV-GFP infection. After exclusion of 

300 cytotoxic siRNA SMARTpools that decreased GFP expression (see above), the remaining 684 targets 

301 were ranked on the basis of the GFP signal in host factor-depleted SARS-CoV-GFP-infected cells 

302 compared to control cells (Fig. 3B). Targets were qualified as antiviral or proviral hits if GFP expression 

303 differed significantly from that in infected control cells transfected with the scrambled siRNA pool (p < 

304 0.05). Knockdown of the majority of the targets (552 proteins) did not significantly alter GFP reporter 

305 gene expression (p > 0.05). However, as is not uncommon in this type of screening experiments and 

306 considering the issue with the origin of the cell line used (see above), we cannot formally exclude that our 
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307 results were influenced to some extent by insufficient knockdown of certain target genes by the siRNA 

308 pools in the library. 

309 Using the criteria outlined above, a total of 90 cellular proteins (19.4% of all targets) were 

310 identified as antiviral factors, since their depletion significantly increased GFP expression, although for 

311 most of them less than two-fold. The ten best antiviral hits are depicted in Fig. 4A and the complete 

312 dataset is provided in Dataset SI. Forty proviral factors were identified and the knockdown of nine of 

313 those reduced GFP expression by more than two-fold (Fig. 4B; for the complete dataset, see Dataset SI). 

314 Although, according to the criteria formulated above (p < 0.05), ANGPT4 (214%; p = 0.0555) and PKR 

315 (210%; p = 0.0884) formally did not qualify as antiviral hits, we have included these proteins in view of 

316 the exceptionally strong stimulation of GFP expression triggered by their knockdown (Fig. 4A). 

317 Furthermore, since its knockdown resulted in an almost 3-fold decrease of the GFP signal (35%; p = 

318 0.0004), DGKE was included as a proviral hit, despite the fact that the viability assay did not rigorously 

319 exclude cytotoxic effects for this siRNA pool (viability 88%, p = 0.0540). 

320 The pro- and antiviral hits identified in the siRNA screen were mapped to cellular pathways using 

321 the IPA software package. Fig. 5 shows the canonical pathways and more general functional categories 

322 (highlighted in color) in which the proviral (green) and antiviral (red) hits were strongly represented (p < 

323 0.05). These pathways included apoptosis, cellular immune response, growth factor signaling, cellular 

324 homeostasis, metabolism of complex lipids, and intracellular and second messenger signaling. 

325 

326 Evaluation of antiviral hits. An unexpectedly large number of antiviral hits was identified in the primary 

327 siRNA screen, although for most of them knock-down resulted in a less than 2-fold increase in SARS- 

328 CoV driven GFP expression. To assess the overall quality of our siRNA screen, and the reliability of the 

329 identification of antiviral hits in particular, we selected a set of strong and weak antiviral hits for further 

330 evaluation, namely PKR, ANGPT4, CLK1 (>2-fold increase in GFP signal), CDK6 (1.8-fold increase), 

331 MAP2K3 (1.6-fold increase), and MAP2K1 (1.2-fold increase). Per target, a deconvoluted set of four 

15 


D 

o 


cfWology JVI Journal of Virology Accepted MdnUSCfTpf Posted OfTlme 


332 individual siRNAs was used for additional knock-down experiments, after which SARS-CoV-GFP 

333 replication was quantified and the knock-down efficiency at the protein level for each target was evaluated 

334 by Western blot analysis (Fig. 6). We checked whether there was strong similarity between the target 

335 sequence of our (human) siRNAs and the corresponding Macaca mulatto sequence, considering the origin 

336 of the 293/ACE2 cells (see above). This was the case for all siRNAs except for siRNA CDK6 #4, which 

337 was therefore excluded from further analysis. Using commercially available antisera against the human 

338 proteins, we were unable to reliably detect endogenous expression of ANGPT4 and CLK1. Therefore, it 

339 remains uncertain whether ANGPT4 and CLK1 are true antiviral or false-positive hits, as we could not 

340 determine knockdown levels and correlate these to effects on virus replication (data not shown). 

341 The mitogen-activated protein kinases (MAPKs) were relatively highly represented among the 

342 antiviral hits (see Fig. 5 and Supplemental Dataset SI) and therefore we included MAP2K1 and MAP2K3 

343 in our secondary evaluation. MAP2K1 was a weak antiviral hit in the primary screen, as its depletion led 

344 to a ~ 1.2-fold increased GFP expression. In validation experiments with individual siRNAs, we also 

345 observed a small, but non-significant increase in SARS-CoV-driven GFP expression (Fig. 6A). Western 

346 blot analysis of siRNA-transfected cells that were infected with wt SARS-CoV revealed poor knockdown 

347 efficiencies and a clear correlation between the level of MAP2K1 and SARS-CoV N protein expression 

348 could not be established (Fig. 6B). The siRNA that gave the best knockdown of MAP2K1 (#2) had no 

349 effect on GFP expression, suggesting that this weak antiviral hit was a false positive in the primary screen. 

350 Knock-down of MAP2K3 resulted in a >1.6-fold increased SARS-CoV-driven GFP expression in the 

351 primary screen (see Dataset SI). In cells transfected with a deconvoluted set of individual siRNAs 

352 targeting the MAP2K3 mRNA we observed a significant increase in GFP expression for 3 out of 4 

353 siRNAs (Fig. 6C). However, in siRNA-transfected cells that were infected with wt SARS-CoV, only 

354 introduction of siRNA #3 led to clearly enhanced SARS-CoV N protein expression (Fig. 6D). The other 

355 siRNAs actually reduced expression of N protein, with an apparent correlation between the remaining 
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356 percentage of MAP2K3 and N protein levels. Based on these results, MAP2K3 could therefore not be 

357 confirmed as an antiviral hit. 

358 We were able to detect expression of cyclin-dependent kinase 6 (CDK6; 1.8-fold increase in GFP 

359 expression in the primary screen) and found that CDK6 siRNAs #1 and #2 reduced protein levels by at 

360 least two thirds (Fig. 6F). Transfection with the same siRNAs significantly enhanced SARS-CoV-GFP 

361 replication (Fig. 6E) and in cells infected with wt SARS-CoV this led to a ~1.5 to 2-fold increase in N 

362 protein levels (Fig. 6F). These results suggest CDK6 to be a bona fide antiviral hit, as its depletion leads to 

363 a moderate but significant and reproducible increase in SARS-CoV replication. 

364 Taken together, the results of our validation experiments suggest that the antiviral hits identified in 

365 the primary screen should be considered with caution, as several of them may have been false-positives, 

366 especially those that had a moderate (but significant) effect in the primary screen. 

367 

368 Validation of PKR as an antiviral factor in SARS-CoV replication. PKR was one of the strongest of 

369 the 90 antiviral hits that were identified in the primary siRNA library screen. In two independent follow- 

370 up experiments with re-ordered PKR-specific siRNA SMARTpools, a more than 2-fold increase in GFP 

371 expression by SARS-CoV-GFP was observed (data not shown), suggesting that PKR is a bona fide 

372 antiviral hit. PKR is a serine/threonine protein kinase that is activated by double-stranded (ds)RNA, a 

373 hallmark of RNA virus infection, and the activated form of PKR blocks translation initiation through elF- 

374 2a phosphorylation (reviewed in (50)). 

375 To further validate the antiviral role of PKR in SARS-CoV replication, a deconvoluted set of four 

376 single PKR-directed siRNAs was used, and transfection of 293/ACE2 cells with three of these siRNAs (# 

377 2, 3, and 4) significantly increased SARS-CoV-driven GFP expression (Fig. 7A; black bars). Cell viability 

378 was slightly reduced after transfection with these PKR-directed siRNAs, in particular using siRNA 2 

379 which caused a 14% reduction in cell viability (Fig. 7A; grey bars). Nevertheless, despite the fact that this 
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380 siRNA adversely affected cell viability, an increase rather than a decrease of SARS-CoV-driven GFP 

381 expression was observed. 

382 Transfection with PKR-specific siRNAs reduced PKR levels in 293/ACE2 cells up to 87% 

383 compared to control cells, depending on the siRNA used (Fig. 7B). To verify that PKR knockdown 

384 increased wt SARS-CoV replication, siRNA-transfected 293/ACE2 cells were infected with wt SARS- 

385 CoV and viral protein expression was analyzed by Western blotting. In line with the effect of PKR siRNA 

386 #2 on 293/ACE2 cell viability (Fig. 7A), cells transfected with this siRNA contained reduced levels of (3- 

387 actin, which was used as loading control (Fig. 7C; lower panel). Transfection with two of the four 

388 individual PKR-directed siRNAs (#2 and 3) clearly increased the expression of SARS-CoV N protein 

389 (Fig. 7C, upper panel), and also led to a ~l-log increase in infectious progeny titers (Fig. 7D). Taken 

390 together, the increases in GFP signal, N expression and infectious progeny titer correlate well with the 

391 magnitude of PKR knockdown, which confirms a strong antiviral role for PKR in SARS-CoV-infected 

392 cells. 

393 

394 Confirmation of a proviral role for protein kinase C iota in SARS-CoV replication. 

395 For the evaluation of the proviral hits of the primary screen (see Fig. 4B) four hits were selected that 

396 caused either a 5-fold reduction (COPB2 and CDK5R2) or a more moderate 2-fold reduction in SARS- 

397 CoV-driven GFP expression (IHPK1 and PRKCi). We were unable to detect endogenous CDK5R2 and 

398 IFIPK1 by Western blot, and therefore could not validate the proviral role of these two host factors (data 

399 not shown). The proviral role of PRKCi and COPB2 could be validated as discussed in the sections below. 

400 The proviral effect of protein kinase C iota (PRKCi) was validated using the chemical inhibitor sodium 

401 aurothiomalate (ATM), which blocks the interaction between PRKCi and other PB1-domain-containing 

402 proteins (51, 52). VeroE6 and 293/ACE2 cells infected with SARS-CoV-GFP were treated with 0.13 to 

403 20pM ATM, starting 2 h prior to infection. Both in 293/ACE2 (Fig. 8A) and in VeroE6 cells (Fig. 8B), 

404 SARS-CoV-mediated GFP expression was efficiently inhibited by ATM in a dose-dependent manner with 
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405 EC 50 values of 0.58 and 1.06 pM, respectively. No cytotoxicity was observed at the ATM concentrations 

406 used (Fig. 8). 

407 

408 COPB2 and proteins of the early secretory pathway are important for SARS-CoV replication. 

409 COPB2 (or P’-COP) was identified as the strongest proviral hit in our screen, as its knockdown resulted in 

410 an 82% decrease of GFP expression (Fig. 4B). The coatomer protein complex, of which COPB2 is a 

411 subunit, contains a total of seven protein subunits (a-, P-, P’-, y-, 8-, 8-, and ^-COP), and drives the 

412 formation of COPI-coated vesicles, which function in retrograde transport in the early secretory pathway 

413 (53). To validate its role as a proviral host factor in SARS-CoV replication, COPB2 was depleted by 

414 transducing 293/ACE2 cells with lentiviruses expressing COPB2 mRNA-specific shRNAs. This reduced 

415 COPB2 levels by ~70%, compared to control cells transduced with a lentivirus expressing a scrambled 

416 shRNA (Fig. 9A), a reduction that did not affect cell viability (Fig. 9B). Subsequent infection of COPB2- 

417 depleted cells with SARS-CoV-GFP resulted in a strong decrease of N protein and GFP expression (Fig. 

418 9C; left panels). To exclude that the observed effect was an artifact caused by the use of the GFP reporter 

419 virus, we also analyzed viral protein expression and virus yield in COPB2-depleted cells infected with wt 

420 SARS-CoV. As for SARS-CoV-GFP, a clear reduction in N protein expression was then also observed in 

421 COPB2-depleted cells, compared to cells transduced with a lentivirus expressing a scrambled shRNA 

422 (Fig. 9C; right panels). Titration of culture supernatants from SARS-CoV-GFP-infected cells and wt 

423 SARS-CoV-infected cells revealed a 2- to 3-log reduction for both viruses upon COPB2 depletion (Fig. 

424 9D). 

425 To further substantiate the importance of COPI-coated vesicles for SARS-CoV replication, 

426 another component of the coatomer protein complex, subunit (31 (COPB1) was depleted by transfection of 

427 293/ACE2 cells with a COPB1 mRNA-specific siRNA SMARTpool. Depletion of COPB1 resulted in a 

428 83% reduction of SARS-CoV-driven GFP expression (Fig. 9E). The formation of COPI-coated vesicles is 

429 mediated through activation of ADP-ribosylation factor 1 (Arfl) by Golgi-specific brefeldin A-resistance 
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430 guanine nucleotide exchange factor 1 (GBF1) (54). Therefore, we also analyzed the importance of GBF1. 

431 GFP reporter gene expression by SARS-CoV-GFP was reduced by 89% in 293/ACE2 cells that had been 

432 depleted for GBF1 (Fig. 9E). GBF1 and COPB1 depletion had no significant effect on cell viability (Fig. 

433 9E). Taken together, these data suggest that COPB2 and COPI-coated vesicles play an essential role in 

434 SARS-CoV replication. 

435 

436 
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437 Discussion 

438 In the past decade, functional genomics studies have systematically identified host factors that can 

439 influence the replication of diverse +RNA viruses (3, 4, 8-10, 39, 55). Here we describe a human kinome- 

440 wide siRNA screen that aimed to identify factors influencing the entry and replication of SARS-CoV. To 

441 our knowledge, this is the first report on a systematic functional genomics study of this kind for any 

442 coronavirus. As kinases are key regulators of many cellular processes, the pro- and antiviral factors 

443 identified in this study should pinpoint cellular pathways that are important for SARS-CoV replication. 

444 After we had completed our screen we unfortunately discovered that the 293/ACE2 cells used 

445 were not of human origin, but must have derived from a non-human primate, probably an Old-world 

446 monkey closely related to Mcicaca mulatto. This may have increased the number of false-negative hits, 

447 due to mismatches between siRNAs designed to target human genes and the sequence of the homologous 

448 monkey mRNAs. However, the human and Macaco mulatto genome are 94% identical and even a 

449 nucleotide mismatch in an siRNA would not automatically render it inactive, as it might still silence gene 

450 expression by blocking translation of the mRNA (56). We therefore concluded that, despite this post- 

451 screening complication, we should still be able to identify host factors that are relevant for coronavirus 

452 infection. Moreover, it is important to stress that the use of a non-human primate cell line should not have 

453 increased the number of false-positive hits. The cell line was highly susceptible to both SARS-CoV and 

454 EAV infection and could be efficiently transfected with siRNAs. This allowed us to perform siRNA 

455 screens for host factors involved in the replication of these two distantly-related ni do viruses, and to 

456 directly compare hits (K.F. Wannee, A.H. de Wilde et al., manuscript in preparation). 

457 A recombinant SARS-CoV-GFP reporter virus, in which ORF7a and ORF7b were replaced by the 

458 GFP gene, was used in our screen in order to conveniently quantify the effect of gene knock-down on 

459 virus replication. The SARS-CoV ORF7a protein is known to interact with the structural envelope (E), 

460 membrane (M), and spike (S) proteins (57), and some studies suggest that it may be involved in specific 

461 virus-host interactions (35, 58, 59). However, the replication efficiency in cell culture of SARS-CoV 

21 


D 

o 


cfWology JVI Journal of Virology Accepted MdnUSCfTpf Posted OfTlme 


462 mutants lacking both ORF7a and ORF7b is unchanged (60, 61), nor were differences in replication 

463 kinetics, morbidity, and mortality observed in a hamster infection model (61). Although the deletion of the 

464 two accessory protein genes likely has affected the results of our primary screen only marginally, a wt 

465 SARS-CoV isolate was used in several of the validation experiments to rule out artifacts caused by the 

466 lack of expression of ORF7a and 7b. In all cases tested, we did not find major differences between wild- 

467 type virus and the deletion mutant. 

468 For SARS-CoV, screening of the kinome-directed library of 779 siRNA SMARTpools resulted in 

469 the identification of 90 antiviral and 40 proviral proteins. Canonical cellular processes and pathways in 

470 which these factors were represented strongly include inositol phosphate metabolism, signaling by Rho 

471 family GTPases, and SAPK/JNK signaling (Fig. 5). Many of the hits could also be mapped to the 

472 interleukin (IL)-2, -6, -8, and IL-17 signaling pathways, which have previously been implicated in 

473 controlling coronavirus infection and coronavirus-induced inflammation (reviewed in (2)). For example, 

474 the SARS-CoV spike (S) protein was shown to induce the expression of the pro-inflammatory cytokine 

475 IL-8 (62), and IL-6 and IL-8 levels were elevated in the serum of SARS-CoV-infected patients (62, 63). 

476 Furthermore, mouse hepatitis virus (MHV) and infectious bronchitis virus (IBV) infections were reported 

477 to upregulate the synthesis of these same cytokines (64, 65). Although our siRNA library screen did not 

478 target interleukins directly, the identification of (kinase-regulated) interleukin signaling pathways is in line 

479 with these earlier studies, and emphasizes their importance in SARS-CoV infection. 

480 A list of host proteins involved in SARS-CoV infection is not more than a good starting point for 

481 follow-up studies into the role of individual host protein or pathways in the vims replication cycle. 

482 Previous studies on other viruses, e.g. HIV-1, showed a very limited overlap between hits from 

483 independent siRNA screens performed in different laboratories (66), highlighting the importance of 

484 validation and follow-up studies. To judge the overall quality of our siRNA screen, six hits with variable 

485 impact on SARS-CoV-driven GFP expression were chosen for validation. Four of these, PKR, CDK6, 

486 COPB2, and PRKQ (Figs. 6-9) could be confirmed. The weak antiviral hits MAP2K1 and MAP2K3 
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487 could not be confirmed, as in follow-up experiments knock-down could not be achieved or did not 

488 convincingly affect SARS-CoV replication, respectively. Interestingly, the diacylglycerol kinase was 

489 highly represented as hits in the primary screen (6 of 8 targets, Supplemental Dataset SI). Although not 

490 included in follow-up experiments presented here, a related study that aimed to identify host factors with a 

491 general effect in nidovirus infection (Fig. 4 and K.F. Wannee, A.H. de Wilde et al., manuscript in 

492 preparation) confirmed that one of these hits, diacylglycerol kinase epsilon (DGKE), plays a role in the 

493 SARS-CoV replication cycle. 

494 MAP2K3, a kinase that acts in the p38 MAPK module, was a moderate hit in our primary screen. 

495 This MAP kinase signaling pathway is involved in multiple processes like regulation of inflammatory 

496 responses, cell proliferation, and cell cycle progression (reviewed in (67). This pathway has been 

497 implicated in the replication of other coronaviruses, but its exact role is still not fully understood. 

498 Activation of p38 MAPK promotes MHV replication (65) and chemical inhibition of the p38 MAPK 

499 pathway restricts HCoV-229E replication (68). Overexpression of the SARS-CoV ORF3a (69) and ORF7a 

500 (35) proteins activates the p38 MAPK signaling pathway, but the role of this pathway in SARS-CoV- 

501 infected cells remains unclear. Our screen identified several proteins from MAPK signaling cascades, but 

502 our validation studies suggested that MAP2K3 was a false positive hit. This is supported by the fact that 

503 the inhibitor SB203580 had no effect on SARS-CoV replication in cell culture (data not shown). This 

504 compound was previously shown to block HCoV-229E infection in cell culture (68) and to increase the 

505 survival of SARS-CoV-infected mice through reducing the SARS-CoV-induced inflammatory response 

506 (70). 

507 CDK6, a kinase involved in cell cycle progression from G1 to S phase (71), was confirmed as an 

508 antiviral hit. Depletion of CDK6 results in G1 phase cell cycle arrest. In addition CDK6 is also involved in 

509 NF-kB signaling and co-regulation of inflammatory genes by binding and activation of the p65 subunit of 

510 NF-kB (72, 73). Consequently, besides the effect on the cell cycle, depletion of CDK6 might also reduce 

511 the inflammatory response against virus infection. A recent study by DeDiego et al. highlighted the 
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512 relevance of NF-KB-mediated inflammation in SARS-CoV-infected mice (74). Several laboratories 

513 studied the (antiviral) role of CDK6, and cell cycle progression in general, in coronavirus replication. For 

514 MFIV, it has been shown that upon high-MOI inoculation, the virus induces cell cycle arrest in the G0/G1 

515 phase to promote its replication. In addition, CDK6 is downregulated in MHV-infected 17C11 cells (75). 

516 Similar observations have been made for SARS-CoV, with the N protein limiting cell cycle progression 

517 by reducing CDK4 and CDK6 kinase activity (76). Overexpression of the SARS-CoV ORF7a protein 

518 induced cell cycle arrest in the G0/G1 phase, however this was not associated with inhibition of CDK4 

519 and CDK6 activity (59). In our study, the antiviral role of CDK6 was confirmed and the observed antiviral 

520 effect is in line with previous studies. 

521 As pointed out above, our screen yielded a relatively high proportion of antiviral hits, although 

522 their effect on SARS-CoV replication, while being statistically significant, was generally limited. Based 

523 on our assessment of some of these moderate hits, at least some of them must have been false-positives. 

524 Knockdown of PKR had the strongest effect (~2-fold increase in GFP expression) on SARS-CoV 

525 replication, and this hit could be confirmed independently, as three out of four individual PKR-directed 

526 siRNAs induced a clear increase in SARS-CoV protein expression and virus yield (Fig. 7C-D). PKR is 

527 one of four mammalian kinases that can phosphorylate eIF-2a in response to stress signals (the others 

528 being the PKR-like endoplasmic reticulum kinase (PERK), GCN2, and HRI). Many virus families have 

529 evolved gene products and strategies to counteract or evade the antiviral action of PKR, highlighting the 

530 importance of this kinase in the antiviral defense. Previously, it was found that PKR inhibits the 

531 replication of the coronavirus IBV, as overexpression of a dominant negative kinase-defective PKR 

532 mutant enhanced IBV replication by almost 2-fold. Furthermore, IBV appeared to (weakly) antagonize the 

533 antiviral activity of PKR through two independent mechanisms, including a partial block of PKR 

534 activation (77). Interestingly, MHV-A59 infection in L2 or 1701 cells did not induce PKR activation and 

535 the sensitivity of MHV to IFN treatment appeared to be PKR-independent (33, 78, 79). TGEV protein 7 

536 was shown to counteract PKR activation by binding protein phosphatase 1 to dephosphorylate eIF-2a 
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537 (80), which - in support of our findings - also suggests that PKR is involved in controlling coronavirus 

538 replication. Krahling et al. showed that PKR was activated in SARS-CoV-infected 293/ACE2 cells, but 

539 concluded that PKR knockdown did not significantly affect virus replication, despite the fact that a ~l-log 

540 increase in SARS-CoV titer was observed in their experiments (81). This is in contrast to our PKR knock- 

541 down experiments, which clearly pointed to an antiviral role for PKR (Fig. 8). In our hands, depletion of 

542 PKR significantly increased SARS-CoV-driven GFP expression (Fig. 8A), and also enhanced N protein 

543 expression (Fig. 8C) and the production of infectious progeny (Fig. 8D). This discrepancy cannot be due 

544 to host cell differences, as the same 293/ACE2 cells were used in both studies (81), and might thus be 

545 attributed to differences in the experimental set-up, choice of controls, or normalization and interpretation 

546 of the data. 

547 In line with the findings for PKR, reducing the expression of PERK (or EIF2AK3), one of the 

548 other kinases known to phosphorylate eIF-2a, resulted in an increase of SARS-CoV-GFP reporter gene 

549 expression with 57% (p < 0.01; Dataset SI). The unfolded protein response - i.e. the detection of 

550 misfolded proteins within the ER lumen - activates PERK, which in turn phosphorylates eIF2a, and 

551 ultimately triggers apoptosis. The relatively strong antiviral effect of PERK observed in this study is in 

552 line with previous studies suggesting that the phosphorylation of eIF-2a in SARS-CoV-infected cells is 

553 mediated by PERK activation (81). Our findings support the hypothesis that upon SARS-CoV infection 

554 the unfolded protein response is activated as an antiviral strategy. In fact, countering SARS-CoV infection 

555 may involve multiple cellular responses that induce apoptosis, including the activation of PKR and PERK, 

556 which could also explain the identification of several other hits involved in apoptosis, like those from the 

557 SAPK/JNK pathway. 

558 Among the proviral hits, PRKCi had a relatively moderate effect in the primary screen (Fig. 4B), 

559 but its proviral role was validated using the inhibitor ATM (Fig. 8). Members of the protein kinase C 

560 family are serine/threonine protein kinases and involved in several signalling pathways that regulate e.g. 

561 cell proliferation, cell cycle progression, and cell survival (reviewed in (82)). Interestingly, PRKCi can be 
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562 activated by phospho-inositol lipids involved in microtubule dynamics within the early secretory pathway 

563 (83). PRKCi contains an N-terminal PB1 domain that ensures the signalling specificity (reviewed in (84)) 

564 and ATM affects the interaction of PRKCi with other PB1 domain-containing proteins like Par6, MEK5, 

565 and p62 (51, 52, 84). Therefore, blocking the PRKCi PB1 domain could decrease MEK5 (85) and NF-kB 

566 signalling (via p62), and affect cell polarity (via Par6) (reviewed in (84)). In addition, PRKCi plays an 

567 essential role in COP-I vesicle formation, since the GTPase Rab2 binds PRKCi and ultimately promotes 

568 recruitment of (3-COP to pre-Golgi membrane structures for the formation of early secretory vesicles (83, 

569 86, 87). As discussed below, COPB2 and the early secretory pathway play a crucial role in SARS-CoV 

570 replication, and in this manner PRKCi may affect SARS-CoV replication as well. Tisdale el al. have 

571 shown that PRKCi kinase activity is essential for the generation of retrograde-transport vesicles (87). 

572 However, the role of the PRKCi PB1 domain, the main target of the drug ATM, in COP-I vesicle 

573 formation was not directly investigated in this study. Nonetheless, this hypothesis is substantiated by the 

574 ATM concentration that blocked SARS-CoV replication. Our EC 5 o values are similar to the reported IC 5 o 

575 (~1 pM) for inhibition of the PRKCi PB-mediated interactions (51), while the IC 5 o for inhibition of the 

576 kinase activity was ~100-fold higher (88). The exact mechanism by which PRKCi influences the SARS- 

577 CoV replicative cycle remains an interesting topic for future research, and it could even be an interesting 

578 target for the development of host-directed antiviral therapy for pathogenic coronaviruses. 

579 Coronavirus replication is associated with a cytoplasmic reticulovesicular network of modified 

580 ER, including double-membrane vesicles and convoluted membranes (28). Despite the in-depth 

581 characterization of their ultrastructure, the biogenesis of these membrane structures and the cellular factors 

582 involved have remained largely uncharacterized. For example, the membrane source of these coronavirus- 

583 induced replication structures is still controversial, with advanced EM analyses showing the RVN to be 

584 derived from and continuous with the ER (28, 44, 89) and other studies implicating the autophagy 

585 pathway (90) or EDEMosomes (91) as the primary membrane donor. Our earlier work already suggested 

586 that the integrity of the early secretory pathway is important for efficient SARS-CoV replication, as 
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587 brefeldin A treatment of SARS-CoV-infected cells significantly reduced replication as well as the 

588 accumulation of virus-induced membrane structures (89). In line with these findings, COPI-coated 

589 vesicles were implicated in the biogenesis of MHV replication structures (92, 93) and SARS-CoV nsp3 

590 was shown to interact with three COPI subunits (94). In none of these previous SARS-CoV and MHV 

591 studies a complete block of virus replication could be achieved, neither by reducing COPI vesicle 

592 formation by depletion of one of the coatomer subunits, nor by treatment with brefeldin A. These results 

593 may in part be explained by incomplete knockdown or the presence of residual COPI vesicles (complete 

594 knockdown is probably not possible due to its detrimental effect on intracellular trafficking and cell 

595 viability). Although our present study clearly demonstrates the importance of COPI-vesicles in SARS- 

596 CoV replication (Fig. 9), their role in the formation or function of the SARS-CoV-induced RVN remains 

597 elusive. The importance of COPI-coated vesicles is further supported by their essential role in the 

598 replication of many other RNA viruses, such as poliovirus (95, 96) and other picomaviruses (40, 97-100). 

599 In conclusion, our kinome-wide siRNA screen has identified several cellular proteins and 

600 pathways that influence SARS-CoV replication and possibly coronavirus infections in general. Our data 

601 thus provide a starting point for further validation and in-depth mechanistic studies which should enhance 

602 our understanding of the complex interplay between coronaviruses and their host. 

603 
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613 Figure legends: 

614 

615 Fig. 1. Viability and susceptibility to SARS-CoV infection of siRNA-transfected 293/ACE2 cells. (A) 

616 293/ACE2 cells were infected with SARS-CoV-GFP (MOI 10) and at 24, 28, and 30 h p.i. cells were 

617 fixed and GFP fluorescence was measured. (B) 293/ACE2 cells were transfected with siRNAs targeting 

618 GAPDH rnRNA or a scrambled control siRNA (Scr). At 48 h p.t., cells were infected with SARS-CoV- 

619 GFP (MOI 10) and 24 h later cells were fixed and GFP expression was measured (black bars). Cell 

620 viability (dark grey bars) was analyzed at 48 h after siRNA transfection and knockdown of GADPH 

621 expression was monitored by measuring enzymatic activity (light grey bars). All values were normalized 

622 to those obtained with non-transfected control cells (100%). 

623 

624 Fig. 2. Design of siRNA library screening procedure. See text for details. 

625 

626 Fig. 3. Results of the siRNA screens for host factors influencing SARS-CoV replication. (A) Viability 

627 assays were done at 48 h p.t. and data were normalized to the measurements for control cells transfected 

628 with scrambled siRNA (100%). Data were binned into 4 viability categories as indicated below the x-as 

629 and the number in each bar is the absolute number of siRNA targets within that category. The fraction of 

630 the total (779) number of targets in each category is indicated above each bar. For each siRNA pool in the 

631 library, the viability data are the average of nine measurements, resulting from three independent library 

632 screens. (B) The plot shows the distribution of the log 2 -transformed values of GFP reporter gene 

633 expression by SARS-CoV-GFP in siRNA-transfected cells, normalized to the GFP signal of infected 

634 control cells that were transfected with scrambled siRNA. Targets were ranked based on the magnitude of 

635 the effect of their knockdown on SARS-CoV replication. Targets were considered to have a robust 

636 antiviral effect when their knockdown increased reporter gene expression to at least 150% (red area above 

637 x-axis). Proviral hits, whose knockdown induced an at least 2-fold reduction in GFP expression, are 
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638 depicted in the green area below the x-axis. Proviral targets whose knockdown reduced cell viability to 

639 below 85% were excluded (see main text), leaving a total of 684 targets included in the final analysis. The 

640 plot represents the average of three library screens (each done in triplicate). 

641 

642 Fig. 4. Heat-maps of pro- and antiviral hits identified in this study. (A) List of the ten most prominent 

643 antiviral (A) and proviral hits (B). For each target, the p-value, accession number, and gene name are 

644 shown. Each data point represents the result of a single library screen and is the average of the 3 replicates 

645 that were done in each screen. The full hit-lists can be found in Dataset SI. 

646 

647 Fig. 5. Cellular pathways influencing SARS-CoV-GFP replication. Graphical representation of the 

648 canonical pathways (white ellipses) identified in the siRNA library screen for cellular factors affecting 

649 SARS-CoV replication. All proviral (green) and antiviral hits (red) (for a complete list, see Dataset S1) for 

650 which depletion significantly altered SARS-CoV-GFP replication were used to identify cellular pathways 

651 by IPA in which the hits were clearly overrepresented, and only those pathways and the hits represented in 

652 them are shown here. The hits are represented by nodes with lines linking them to one or more canonical 

653 pathways. The color intensify of the nodes indicates the strength of the pro- or antiviral effect (log 2 -ratio of 

654 GFP expression normalized to infected control cells), e.g. factors with a stronger antiviral effect have a 

655 more intense red color. The identified canonical pathways were clustered into more general categories that 

656 are indicated by text boxes in the colored background shading. 

657 

658 Fig. 6. Evaluation of the antiviral hits CDK6, MAP2K1, and MAP2K3. 293/ACE2 cells were 

659 transfected with four individual siRNAs targeting, MAP2K1 (A, B) or MAP2K3 (C, D), or three CDK6- 

660 specific siRNAs (E, F). A non-targeting scrambled siRNA was used as a control. At 48 h p.t. cells were 

661 infected with SARS-CoV-GFP at an MOI of 10 (A, C, E), fixed 24 h later, and GFP fluorescence (black 

662 bars) was quantified and normalized to the value measured in infected, scrambled siRNA-transfected cells 
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663 (100%). The effect of siRNA transfection on cell viability was analyzed in parallel (grey bars) and values 

664 were normalized to those of scrambled siRNA-transfected control cells (100%). Each experiment was 

665 repeated at least three times (average ± SD). In parallel, siRNA-transfected cells were infected with wt 

666 SARS-CoV (MOI 5), and at 8 h p.i., SARS-CoV N expression was monitored by Western blotting (B, D, 

667 F). TfR was used as a loading control. Knockdown levels of the host proteins were analyzed by Western 

668 blot (B, D, F). The amount of SARS-CoV N protein and remaining quantity of host protein compared to 

669 scrambled siRNA-transfected cells (100%) is shown below each lane. All experiments were repeated at 

670 least twice. 

671 

672 Fig. 7. Validation of PKR as an antiviral factor in SARS-CoV replication. 293/ACE2 cells were 

673 transfected with four individual siRNAs targeting PKR or a scrambled control siRNA. (A) At 48 h p.t. 

674 cells were infected with SARS-CoV-GFP (MOI 10), fixed 24 later, and GFP fluorescence (black bars) was 

675 quantified and normalized to the value measured in infected, scrambled siRNA-transfected cells (100%). 

676 The effect of siRNA transfection on cell viability was analyzed in parallel (grey bars) and values were 

677 normalized to those of scrambled siRNA-transfected control cells (100%). Average ± SD is given (***; p- 

678 value < 0.001). (B) Knockdown of PKR expression at 48 h p.t. was monitored by Western blotting and the 

679 percentage of PKR remaining compared to scrambled siRNA-transfected cells is shown below each lane. 

680 TfR was used as loading control. (C) Cells transfected with PKR-specific siRNAs and control cells were 

681 infected with SARS-CoV (MOI 0.01) and 24 h later these cells were lysed to assess SARS-CoV N levels 

682 by Western blotting (shown below the panels as percentage of control), using (3-actin as loading control. 

683 (D) Vims titers in the 24-h p.i. culture supernatants of wt SARS-CoV-infected cells (MOI of 0.01) 

684 transfected with PKR-specific or scrambled siRNA. All experiments were repeated at least twice. 

685 

686 Fig. 8. Validation of PRKCi as proviral host factor. (A) 293/ACE2 or (B) VeroE6 cells in 96-well 

687 plates were infected with SARS-CoV-GFP (MOI 10). Treatment with 0.13-20 pM sodium aurothiomalate 
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688 (ATM) was started 2 h prior to infection and cells were fixed at 18 h p.i. (VeroE6) or 24 h p.i. 

689 (293/ACE2). GFP reporter gene expression (black bars) was measured and normalized to the signal in 

690 untreated control cells (100 %). The grey lines show the effect of ATM on the viability of mock-infected 

691 cells, normalized to the viability of solvent-treated control cells. Graphs show the results (average and SD) 

692 of a representative experiment performed in quadruplicate. Both experiments were repeated at least twice. 

693 

694 Fig. 9. Proteins of the early secretory pathway are important for SARS-CoV replication. (A) 

695 293/ACE2 cells were transduced with lentiviruses expressing a COPB2 mRNA-specific or a scrambled 

696 shRNA. Knockdown of COPB2 expression at 48 h p.t. was monitored by Western blotting with a COPB2- 

697 specific antiserum and cyclophilin B (CypB) was used as loading control. (B) Viability of COPB2- 

698 depleted 293/ACE2 cells was analyzed at 48 h after transduction (% of control cells transduced with 

699 lentiviruses expressing a scrambled shRNA). (C, D) COPB2-depleted and control cells were infected with 

700 either SARS-CoV-GFP or wt SARS-CoV (MOI of 0.01). (C) SARS-CoV protein expression at 32 h p.i. 

701 (SARS-CoV-GFP) or 24 h p.i. (wt SARS-CoV) was analyzed by Western blotting with N-specific and 

702 GFP-specific antisera, using the TfR protein as loading control. SARS-CoV N protein expression was 

703 quantified and normalized to that inscrambled siRNA-transfected cells (100%) as indicated under each 

704 lane. (D) SARS-CoV-GFP (black bars) and wt SARS-CoV (grey bars) progeny titers in the culture 

705 supernatants of control or COPB2-depleted cells at 32 h p.i. (SARS-CoV-GFP) or 24 h p.i. (wt SARS- 

706 CoV). (E) Normalized GFP expression by SARS-CoV-GFP in 293/ACE2 cells transfected with siRNA 

707 SMARTpools targeting COPB1, GBF1, or a scrambled control siRNA. Cells were infected 48 h p.t. at an 

708 MOI of 10 and 24 h later GFP fluorescence was quantified and normalized to that in infected cells 

709 transfected with a scrambled siRNA. GFP fluorescence data is the average of three independent 

710 experiments. 
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Figure 7 
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